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ABSTRACT 

We present cosmological constraints from 2D weak gravitational lensing by the large- 
scale structure in the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS) which 
spans 154 square degrees in five optical bands. Using accurate photometric redshifts and 
measured shapes for 4.2 million galaxies between redshifts of 0.2 and 1.3, we compute the 
2D cosmic shear correlation function over angular scales ranging between 0.8 and 350 ar- 
cmin. Using non-linear models of the dark-matter power spectrum, we constrain cosmolog- 
ical parameters by exploring the parameter space with Population Monte Carlo sampling. 
The best constraints from lensing alone are obtained for the small-scale density-fluctuations 
amplitude erg scaled with the total matter density O m . For a flat ACDM model we obtain 
cr 8 (fi m /0.27) - 6 = 0.79 ± 0.03. 

We combine the CFHTLenS data with WMAP7, BOSS and an HST distance-ladder prior 
on the Hubble constant to get joint constraints. For a flat ACDM model, we find f2 m = 
0.283 ± 0.010 and <j 8 = 0.813 ± 0.014. In the case of a curved wCDM universe, we obtain 
n m = 0.27 ± 0.03, a 8 = 0.83 ± 0.04, w Q = -1.10 ± 0.15 and n K = 0.006±g;gg£. 

We calculate the Bayesian evidence to compare flat and curved ACDM and dark-energy 
CDM models. From the combination of all four probes, we find models with curvature to 
be at moderately disfavoured with respect to the flat case. A simple dark-energy model is 
indistinguishable from ACDM. Our results therefore do not necessitate any deviations from 
the standard cosmological model. 

Key words: cosmological parameters - methods: statistical 
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1 INTRODUCTION 

Weak gravitational lensing is considered to be one of the most pow- 
erful tools of cosmology. Its ability to measure both the geometry 
of the Universe and the growth of structure offers great potential to 
obtain constraints on dark energy and modified gravity. Moreover, 
to first order, weak lensing does not rely on the relation between 
galaxies and dark matter (bias), and is therefore a key probe of the 
dark Universe. 

Cosmic shear, denotes the distortion of images of distant 
galaxies due to the continuous deflection of light bundles prop- 
agating through the inhomogeneous Universe. The induced cor- 
relations between shapes of galaxies are directly related to the 
statistical properties of the total (dark + luminous) large-scale 
matter distribution. With an estimate of the redshift distribu- 
tion of the lensed galaxies, theoretical predictions of weak- 
lensing observables can be tested to obtain constraints on cos- 
mological parameters and models. Re cent reviews which also 
summ ar ize past observational resu l ts arelBartelmann & SchneideJ 
d200lh: IVan Waerbeke & Mellied (120031); iMunshi et al l d2008l) : 
iHoekstra & Jainl d2008l) : lBartelmannl d2010l) . 

The Canada-France Hawaii Telescope Legacy SurvejQ 
(CFHTLS) is a large imaging survey, offering a unique combina- 
tion of depth (zab Si 24.5 at 5a point source limiting magnitude) 
and area (154 square degrees). It is the largest survey volume over 
which cosmic shear has ever been measured. This paper presents 
the first cosmological analysis of the complete CFHT Legacy Sur- 
vey with weak gravitational lensing. We measure the seond-order 
cosmic-shear functions from the CFHTLenS (CFHT Lensing Sur- 
ve)0), which comprises the final CFHTLS data. Earlier analyses 
of the CFHTLS used the fir st data release (T0001) with 4 square 
degrees of the Deep survey dSemboloni et al.|[2005h and 22 s quare 
degrees of the Wide part of the survey dHoekstra et al.ll200oT) . fol- 
lowed by the third data relea s e (T0003) com prising 55 square de- 
grees dBeniamin et"al] 120071 : iFu et all 120081, hereafter F08). The 
T0003 lensing data was subsequently employed in further studies 
dDore et aljliooil iTereno et al.ll2009t iKilbinger et ai1l2009l) . Only 
CFHTLS-Wide i'-band data were used for those lensing analyses, 
and the redshift distribution w as inferred from th e photometric red- 
shifts from the Deep survey dllbert et al.ll2006l) . Photometric red- 
shifts in the Wide were ob tained subsequently with the T0004 re- 
lease dCoupon et alj|2009h . The current series of papers uses the 
final CFHTLenS data release of 154 square degrees in the five op- 
tical bands u* , g ,r' ,i' , z' . The analysis improved significantly in 
several ways: 

• The CFHTLS data have been reanalysed with a new pipeline 
dErbenetal .1120091 l2013h . 

• Photometric redshifts have been obtained for each individ- 
ual galaxy in the lensing catalogue from PSF-homogenized images 
dHildebrandt et al.ll2012l) . The accuracy h as been verified in det ail 
by an angular cross-correlation technique dBeniamin et al 1120131) . 

• PSF modeling and galaxy shape measurement have been 
performed with the forward model-fitting method lensfit, which 
has been t horoughly tested on simulations and improved for 
CFHTLenS jMiller et al .1120131) . 

• Systematics tests have be en performed in a blin d way, to yield 
unbiased cosmological results dHevmans et al.120121) . 

• The cosmology-dependent covariance matrix is obtained by 

1 http : //www . cf ht . hawaii . edu /Science /CFHTLS 

2 http://www.cfhtlens.org 



a mixture of numerical simulations on small scales and analytical 
predictions on large scales. 

The reliability and accuracy of our photometric redshifts al- 
lows for three-dimensional weak lensing analyses. The measure- 
ment of the lensing correlations for different redshift combinati ons 
allow s us to obtain information on the growth of structure (e.g. iHul 
1 19991) . and has a great po tential to constrain dark energy a nd mod- 
ified gravity models (e.g. lAlbrecht et al .1120061 : IUzanll2010l) . In sev- 
eral companion papers, we perform 3D cosmic shear analyses by 
splitting up galaxies into redshift bins using correlation function 
methods presented here (lensing tomograph y: ISimpson et al.l2013l : 
iHevmans et al.ll2013l : lBeniamin et ai]|2013h . 

In this paper, we perform a 2D lensing analysis using a single 
redshift distribution. Despite the fact that the redshift information 
is not used in an optimal way, our analysis has several advantages. 
First, it yields the highest signal-to-noise ratio (S/N) for a sin- 
gle measurement. This is particularly important on large angular 
scales, where the S/N is too low to be used for tomography. These 
large scales probe the linear regime, where non-linear and bary- 
onic effects do not play a role, a nd one can therefore ob tain very 
robust constraints on cosmology dSemboloni et aljfeoi ll) . Second, 
we can include low-reds hift galaxies without h aving to consider 
intrinsic alignments (IA; iHirata & Seliakj|2004l) . For a broad red- 
shift distribution, IA is expected to be a sub-dominant contribution 
to the cosmological shear-shear correlation with an ex pected bias 
for as which is well within the statistical uncertainty dKirk et al] 
l20ld : iMandelbaum et alj|201ll : Ijoachimi et alj|2oTH) . see also a 
joint lensing a nd IA tomography an alysis over the full available 
redshift range dHevmans et al .1 KO 1 3l) . Therefore, despite the fact 
that a 2D lensing is more limited than tomography, it is less noisy 
and more immune to primary astrophysical systematics. It is there- 
fore a necessary basic step and puts any further cosmological ex- 
ploitation of CFHTLenS using more advanced tomographic or full 
3D lensing techniques on solid grounds. Such analyses are pre- 
sented in the CFHTLenS companion papers. 

This paper is organized as follows. Section [2] provides the 
expressions for the second-order lensing observables used in this 
analysis, both obtained from theoretical predictions and estimated 
from data. The measured shear functions and covariances are pre- 
sented in Sect. [3] Cosmological models and sampling methods are 
introduced in Sect. [4] The results on cosmological parameters and 
models are presented in Sect. [5] followed by consistency tests in 
Sect. [6] The paper is concluded with a discussion in Sect. [7] 

2 WEAK COSMOLOGICAL LENSING 

In this section the main relations between second-order weak 
lensing observables and cosmological qu antities are given. See 
iBartelmannl d20ld) : lHoekstra & Jainl d2008l) for recent reviews. 

2.1 Theoretical background 

Weak lensing by the large-scale structure measures the conver- 
gence power spectrum P K , which can be related to the total mat- 
ter power spe ctrum P$ via a projection using Limber's equation 
dKaiseJl992T) : 

The projection integral is carried out over comoving distances x> 
from the observer out to the limiting distance xiim of the survey. 
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The lens efficiency G is given by 



G(X) = 5 ( — 



/ ,ta dxW)^, (2) 
o(x) J x Ik [X ) 



where Ho is the Hubble constant, c the speed of light, f2 m the total 
matter density, and a(x) the scale factor at comoving distance \- 
The comoving angular distance is denoted with Jk and depends on 
the curvature K of the Universe; 



' K~^ 2 sin (K'^w) 



for A > 

/k (w) = ^ w for A" = 

{-K)- 1/2 sinh ((-A) 1/2 ™) for A' < . 

(3) 

The 3D power spectrum is evaluated at the wave number k — 
l/fK{x)> where I denotes the projected 2D wave mode. The func- 
tion p represents the weighted distribution of source galaxies. 



2.2 Flavours of real-space second-order functions 

From an observational point of view, the most direct measurement 
of weak cosmological lensing is in real space, by using the weak 
gravitational shear signal as derived from galaxy ellipticity mea- 
surements. The two-point shear correlation functions (2PCFs) 
and are estimated in an unbia sed way by averaging over pairs 
of galaxies dSchneider et al.l2002h . 



Eij WiWj [et(tfQ et(flj) ± e x (flQ ex Cgjjj 



(4) 



The sum is performed over all galaxy pairs (ij) with angular dis- 
tance \&i — within some bin around ■&. With et and £ x we 
denote the tangential and cross-component of the galaxy ellipticity, 
respectively. The wei ghts Wj are obtain ed from the /e/7.?fit shape 
measurement pipeline jMilleretal.120131) . The 2PCFs are the Han- 
kel transforms of the convergence power spectrum P K or, more pre- 
cisely, of linear combinations of the E- and B-mode spectra, Pe and 
Pb > respectively. Namely, 



1 

2^ 



All [Pe(£)±Pb(£)}Jo A (M), (5) 



where Jo and J4 are the first-kind Bessel functions of order and 
4, and correspond to the components and , respectively. 

It is desirable to obtain observables which only depend on the 
E-mode and B-mode, respectively. Weak gravitational lensing, to 
first order, only gives rise to an E-mode power spectrum and there- 
fore, a non-detection of the B-mode is an important sanity check of 
the data. To this end, we calculate the following second-order shear 
quantities which can be derived fro m the correlation funct ions: 
The aperture-mass dispe rsion (Mf D ) dSchneider et al.|[l998h . the 
shea r top-hat rms (h \ 2 ) dKaiseJll992h . the optimized ring statistic 
TZe dFu & Kilbingetold) and COSEBIs (Comp lete Orthogonal 
Sets of E-/B-mode Integrals; ISchneider et alj2010l) . The optimized 
ring statistic was introduced as a general isation of the so-called 
ring statistic dSchneider & Kilbingeg|2007l) . The corresponding fil- 
ter functions have been obtained to maximise the figure-of-merit of 
fi m and erg for a CFHTLS-T0003Tike survey. We use these func- 
tions for CFHTLenS which, despite the larger area, has similar sur- 
vey characteristics. COSEBIs represent yet another generalisation 
and contain all information about the E- and B-mode weak-lensing 
field from the shear correlation function on a finite angular range. 

Being quantities obtained from the 2PCFs by non-invertible 
relations, these derived functions do not co ntain the full info rma- 
tion about the convergence power spectrum dEifler et al.ll2008l ). but 



separate the E- and the B-mode in a more or less pure way, as will 
now be described. The derived second-order functions can be writ- 
ten as integrals over the filtered correlation functions. They can be 
estimated as follows: 

A e ,b = ~ J2 ft Atfi [F+ (*<) e+(tfc) ± P_ (0<) £_ (A)] . (6) 

i 

Here, A$j is the bin width, which can vary with i, for example in 
the case of logarithmic bins. With suitable filter functions F+ and 
F- (App.IaI. the estimator Xe (Xb) is sensitive to the E-mode (B- 
mode) only. The filter functions are defined for the various second- 
order observables in Table [TJand Appendix [A] 

All derived second-order functions are calculated for a family 
of filter functions. For the aperture-mass dispersion, the optimized 
ring statistic, and the top-hat shear root mean square (rms), these 
are given for a continuous parameter 9 which can be interpreted 
as the smoothing scale. Here and in the following we will use the 
notation as the scale for the 2PCFs, and '6' as the smoothing 
scale for derived functions. 

For COSEBIs, the filter functions are a discrete set of func- 
tions. The latter exist in two flavours, Lin-COSEBIs and Log- 
COSEBIs, defined through filter functions F± which are polyno- 
mials on linear and logarithmic angular scales, respectively. Here 
we use Log-COSEBIs, for which many fewer mo des are required 
to capture the same info rmation as Lin-COSEBIs dSchneider et alj 
l20ld : lAsgari et alj|2012l) . See App.|A]for more details. 

All of the above functions can be expressed in terms of the 
convergence power spectrum. The general relation is 



XE ' B = 2^ 



dllP E , B (l)U 2 (l). 



(7) 



The functions F± and U 2 are Hankel- transf orm pairs, their 
relatio n is given by ICrittenden et al.l J2002I) and ISchneider et alj 
d2002l) as 



F±(x)= / dtt 3 0A (xt)U 2 {t). 
Jo 



(8) 



2.2.1 Finite support 



The measured shear correlation function is available only on a finite 
interval [# m i n ; $max]- The upper limit is given by the finite survey 
size. We choose 460 arc minutes, which is roughly the largest scale 
for which a sufficient number of galaxy pairs are available on at 
least three of the four Wide patches[f]. The lower limit comes from 
the fact that for close galaxy pairs, shapes cannot be measured reli- 
ably. For lensfii, galaxies with separations smaller than the postage 
stamp size of 48 pixels ~ 9 arcsec tend to ha ve a shape bias along 
their connecting direction dMiller et alj|2013l) . Since this direction 
is randomly orientated to a very good approximation, we do not 
have to remove those galaxies altogether; setting a minimum an- 
gular separation of i9 m i n = 9 arcsec does avoid the correlation of 
those pairs. Both galaxies from a close pair can be independently 
correlated with other, more distant galaxies, for which the close- 
pair shape bias acts as a second-order effect and can safely be ne- 
glected. 

If the support of the filter functions F+ and P_ exceeds the 
observable range, eq. ((6) leads to biased results, and a pure E- and 
B-mode separation is no longer guaranteed. 



W2 as the smallest patch extends to 400 arc minutes, whereas Wl probes 



scales as large as 685 arc minutes. 
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Table 1. E- and B-mode separating second-order functions. They are estimated as integrals, or sums in the discrete case, over the 2PCFs £± multiplied with 
the filter functions F± with the formal integration range [i9 m j n ; tax], see eq. (6). The argument ■& denotes the integration variable. 



Name 




F± (eq.g) 


U (eq.|7j 


"mm 


^max 


Reference 


Aperture-mass dispersion 




T ± (ti/9)/8 2 (eq.ED 


eq. (A2) 





26 


Schneider et al. (1998) 


Top-hat shear rms 




S±(0/O)/9 2 (eq.|A3) 


eq. <A4l 





CO 


Kaiser (1992) 


Optimized ring statistic 


7e B ,B(f = i? max ) 


T ± (i9)(eq.E5J 


n/a 


u mm 


"max 


Fu & Kilbinger (2010) 


COSEBIs 






n/a 


"mm 


^max 


Schneider et al. (2010) 



This is the case for the aperture-mass dispersion and the top- 
hat shear rms. For the former, only the lower angular limit is prob- 
lematic and causes leakage of the B-mode into the E-mode signal 
on small smoothing scales. On scales larger than 9 mm = 5.5 arc 
min utes however, this leak age is below 1.5 cent for i9 m in = 9 arc- 
sec l lKilbinger et alj[200fj) . We therefore choose # m m = 5.5 ar- 
cmin to be the smallest smoothing scale for (M 2 p){9). Note that 
on scales smaller than $ m in we set the correlation function to zero 
and do not use a theoretical model to extrapolate the data on to this 
range, to avoid a cosmology-dependent bias. 

The B-mode leakage for the top-hat shear rms ( \ 'y 2 | ) is a func- 
tion of both the lower and upper available angular scales. Over our 
range of scales, the predicted leaked B-mode for the WMAP7 cos- 
mology is nearly constant with a value of 5.3 x 10~ 7 . 

The first pure E-/B-mode separating function for which the 
c orresponding filter functions h ave finite support was introduced 
in Schneide r & KilbingeJ d2007h . Following that approach, the op- 
timized ring statistic and COSEBIs were constructed in a similar 
way to not suffer from an E-/B-mode leakage. 

An additional bias arises from the removal of close galaxy 
pairs i n the lensing analysis, as was first reported bv lHartlap et alj 
d201ll) . As previously discussed, lensfil produces a shape bias for 
galaxies separated by less than 9 arcsec, but this bias is random and 
the close pairs are therefore used in the analysis, to be correlated 
with other galaxies at larger distances. That said, for very close 
blended galaxies, where it is non-trivial to determine if there is one 
or two galaxies observed, galaxy shapes cannot even be attempted. 
These blended pairs are therefore not reliably detected and lost 
from our analysis. This causes a potential bias, since these galaxies 
are removed preferentially at low redshift, where galaxy sizes are 
larger, and from high-density regions compared to v oids, because 
alaxi es trace the large-scale structure. From Fig. 3 of lHartlap et alj 
20 111) we infer that the magnitude of this effect is at the per cent 
level on scales larger than 0.8 arcmin for the 2PCFs. 



3 CFHTLENS SHEAR CORRELATION DATA AND 
COVARIANCE 

The CFHTLenS data are described in severa l companion papers; 
for a full summary see iHevmans et alj 1 2012b. Stringen t systemat- 
ics tests have been performed in Hevmans et alj |2012) which flag 
and remove any data in which significant residual systematic s are 
detected. It is this cleaned sample, spanning ~ 75 per cent of the 
total CFHTLenS survey area, that we use in this analysis. This cor- 
responds to 129 out of 171 MegaCam pointings. In this paper, we 
complement those tests by measuring the B-mode up to large scales 
(Sect. 13.51 ). Comparing this to the previous analysis of F08 shows 
the improved quality of the lensing analysis by CFHTLenS. 




sum of pdf(z„) 
best-fit z p 



2.5 



3.5 



Figure 1. The redshift distribution p(zp) histogram, estimated as the sum 
of pdfs for 0.2 < z p < 1.3 (thick red curve). For comparison, the his- 
togram of the best-fitting photo-zs is shown as thin blue curve. Error bars 
correspond to the variance between the four Wide patches. 



3.1 Redshift distribution 

A detailed study of the reliability of our photometric redshifts, the 
contaminations between redshift bins, and the c osmological impli- 
cations is performed in iBeniamin et alj d2013l) . This work shows 
that the true redshift distribution p(z) is well approximated by the 
sum of the probability distribution functions (pdfs) for all galaxies. 
The pdfs are output by BPZ (Bayesian Photometric Redshift Es- 
timation; [BenftegHcXx]) a s a function of photomet ric redshift z p , 
and have been obtained bv lHildebrandt et all J2012h . The resulting 
p(z) is consistent with the contamination between redshift bins as 
estim ated by an angular cross-correlation analysis feeniamin et alj 
120101) . The contamination is relatively low for galaxies selected 
with 0.2 < z p < 1.3, which is confirmed by a galaxy -gal axy- 
lensing redshift scaling analysis in lHevmans et al.1 < !2012h . The re- 
sulting p(z) is shown in Fig.Q] The mean redshift is z — 0.748. In 
contrast, the mean redshift of the best-fitting z p histogram is biased 
low with z = 0.69. 



3.2 Angular correlation functions 

We calculate the two-point shear correlation functions by averag- 
ing over pairs of galaxies, using the tree code ATHEN/0. Galaxies 
are partitioned into nested branches of a tree, forming rectangular 
boxes in right ascension a and declination S. For two branches at 
angular distance d and box sizes di (i — 1,2 along the a and 8 di- 
rection, respectively) if both opening angles uii = di /& are smaller 



http : //www2 . iap . f r/users/kilbinge/athena 
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North pole (5 = n/2) 




Figure 2. Angles and coordinates on a sphere for two galaxies i = 1,2 
located at (a;, <5i). Great circle segments are drawn as bold lines. 



than a threshold angle ui t h, the tree is not followed further down by 
descending into sub-branches. Instead, the weighted average shear 
in each branch is used for the 2PCF estimator (eq. |4). We found 
that a value of wth = 0.03 gives sufficient accuracy compared to 
the brute-force approach. 

We use distances and angles on the sphere to calculate the 
shear correlation functions. For two galaxies i — 1, 2 at right as- 
cension and declination {pti, Si), we calculate the great-circle dis- 
tance $ with 



COS 1? = COs(a2 — Ql) COS <5l COS 82 + SHI <5l 



sin 82 



(9) 



Each galaxy's ellipticity is measured in a local Cartesian coordinate 
system with the x-axis going along the line of constant declination 
and the j/-axis pointing to the North pole. We project this ellipticity 
to the tangential and radial component with respect to the connect- 
ing great circle. For that, we calculate the angle fa between the 
great circle segments 1? and on. Then, the projection or so-called 
course angles are <j>i — tt/2 — fa. With the sine and cosine rules on 
the sphere, we get 



cos <pi — 



sin(a2 —0:1) cos 82 
sin ■d 



cos 82 sin <5i — sin 82 cos <5i cos(q?2 — oei) 

sin 0i= : — , (10) 

sin # 

and corresponding expressions for 4>2 by exchanging indices. 

To estimate the smoothed second-order quantities, we com- 
pute the 2PCFs on 10,000 linear angular bins. This is large enough 
not to cause a significant E-/B-mode leakage due to the approxi- 
mation of the integrals over the correlation functions by the direct 
sum (eq. [6j. We verified this using CFHTLenS numerical simula- 
tions with no B-mode (see next Section) , see also lBecketl d2012h . 
We choose the smallest angular distance between two galaxies to 
be 9 arcsec, corresponding to the first bin centre to be 10.4 arcsec. 
We calculate the 2PCFs (eq.[4} as the weighted mean over the four 
Wide patches, using the number of galaxy pairs as weight for each 
bin. 
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Figure 3. Diagonal of the covariance C++ (top panel) and C — (bottom), 
split up into various terms: Shot-noise D (solid red line), mixed term M 
(dashed green), cosmic variance V (dotted blue line and crosses) and shear 
calibration (see Sect.ET4l. 



3.3 Data covariance 

To model and interpret the observed second-order shear functions, 
we need to estimate the data covariance and its inverse. The cos- 
mic shear covariance C is composed of the shot-noise D, which 
only appears on the diagonal, a cosmic-v ariance contribution V, 
and a mixed term M dSchneider et alj|2002h . The covariance of the 
2PCFs is comprised of four block matrices. The diagonal consists 
of C++ and C which are the auto-correlation covariance matri- 
ces of £+ and respectively. The off-diagonal blocks are C+- 
and C-+ = C+_ which denote the cross-correlation covariance 
between £+ and . 

Since the cosmic shear field is non-Gaussian on small and 
medium angular scales, the cosmic variance involves four-point 
functions. Neglecting those can yield overly optimistic cosmo- 
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Figure 4. Correlation coefficient of the total covariance, shot-noise plus 
mixed plus cosmic variance. The lower right triangle (within the red bound- 
ary) is the Gaussian prediction, the upper left (blue) triangle shows the to- 
tal covariance with non-Gaussian parts from the Clone. The four blocks 

are C+_, C , C++ and C |-, respectively, as indicated in the plot. The 

axes labels are the matrix indices % and j . 



logical constraints dSemboloni et al.l |2007| ; iTakada & Jainl |2009| ; 
lHilbertetalj201lh . 

To ac count for non-Gaussianity, we use iV-body simula- 
tions from iHarnois-Deraps et alj d2012l) . From these simulations, 
a 'Clone' of the CFHTLenS data has been produced with the same 
galaxy redshift distribution, galaxy clustering, masks, and noise 
properties. The cosmological lensing signal is added using ray- 
tracing through the light cones. The Clone cosmology is a flat 
ACDM model with f2 m = 0.279, fi b = 0.046, n s = 0.96, <t 8 = 
0.817 and h = 0.701. The lensing signal for each galaxy is con- 
structed by ray-shooting through the simulated dark-matter dis- 
tribution. Each simulated line of sight spans a field of view of 
3.5 x 3.5 square degrees. We fit close to 4 x 4 MegaCam point- 
ings on each line of sight, which is possible because of overlap- 
ping areas between pointings. A total of 184 independent lines of 
sight are used to calculate the field-to-field covariance matrix. The 
final matrix is scaled with the ratio of the effective areas (includ- 
ing masks) of 0.11 which corresponds to 90 per cent of the area 
of 16 MegaCam pointing that fit into each line of sight, divided 
by 129 MegaCam pointings used in this analysis. We average over 
three different samples of the galaxy redshift probability distribu- 
tion, where galaxy redshifts were drawn from the corresponding 
pdf. 

As shown in the upper pa nel of Fig. [3] the Gaussian pre- 
diction for the cosmic variance dKilbinger & Sch neider 200^) for 
£+ provide a good match to the Clone covariance on intermedi- 
ate scales, 10 arcmin S; $ S; 30 arcmin. On larger scales, up 
to i? < 200 arcmin, the numerical simulations under-pred ict the 
power due to the finite box size (e.g. lPower & Knebdl2006l) . Only 



the last two data points show an increased variance, which is due to 
the finite Clone field geometry. When comparing the Clone mean 
correlation function to a theoretical prediction with cut-off scale 
k — (2-7r/147) hT 1 Mpc, we get a rough agreement between the 
two, indicating that the lack of power is indeed caused by the finite 
box. We draw similar conclusions for the cosmic variance of 
shown in the lower panel of Fig. [3] Further, we verified that a Jack- 
knife estimate of the variance by sub-dividing the CHFLTenS data 
into 129 subfields gives consistent results. 

3.3.1 Grafting the covariance matrix 

We construct the total covariance out to ■& — 350 arcmin 
by grafting the Clone covariance V c i to the analytical Gaus- 
sian prediction. For th e latte r, we use the method developed in 
iKilbinger & Schneiderl d2004l) . which takes into account the dis- 
crete nature of the galaxy distribution and the field geometry. First, 
we add the Clone covariance V c i to the Gaussian cosmic covariance 
term Vg . The combined cosmic covariance is 



Vij = gijV \,ij + (1 - gij)VG,ij, 



(11) 



where the modulation function g^ alleviates discontinuities in the 
combined matrix. We choose to be a bi-level step function, with 
<?s,s = 1/2; gij = 1 if both indices i,j are smaller than the step 
index s; and g^ = if at least one of the indices i or j is larger than 
or equal to s. The step index s is chosen such that # s is the scale 
closest to 30 arcmin. Equation Jilt is applied to all covariances 
between the two shear correlation functions, i.e. V++,V+_ and 
V . 

The Clone covariance also cont ains an additiona l variance 
term, which was discovered recently dSatoetal.ll2009l) . This so- 
called halo sample variance (HSV) stems from density fluctua- 
tions on scales larger than the (finite) survey size that are corre- 
lated with fluctuations on smaller scales. For example, the num- 
ber of halos in the survey depend on the large-scale modes out- 
side the survey, since halos are clustered and do not just follow a 
Poisson distribution. This introduces an extra variance to the mea- 
sured power spectrum. The halo sample variance is proportional t o 
the rms density fluctuations at the survey scale dSato et al.ll2009l) . 
Since our simulated light-cones are cut-outs from larger boxes of 
size L = 147 Mpc/h (L = 231 Mpc/h) at redshift below (above) 
unity, they do contain Fourier scales outside the survey volume and 
their coupling to smaller scales. The halo sample variance is im- 
portant on small scales, where our cosmic va r iance is dominated 
by the Clone covariance. Following lSato et all d2009h . we estimate 
the halo sample variance to dominate the CFHTLenS total covari- 
ance at £ « 2 x 10 3 , corresponding to 5 arc minutes which is the 
Clone covariance regime. 

The missing large-scale Fourier modes in the simulation box 
cause the HSV to be underestimated. A further underestima- 
tion comes from the rescaling of the Clone lines of sight to the 
CFHTLenS area since, in contrast to the other covariance terms, 
the HSV term de creases less str ongly than the inve rse survey area 
dSato et alj|2009h . According to iKavo et al] d2012l) . when naively 
rescaling from a 25 deg 2 -survey to 1500 deg 2 , the signal-to-noise 
ratio is too optimistic by not more than 10 per cent. For a re-scaling 
to the smaller CFHTLenS area, this bias is expected to be much 
less. 



3.3.2 Cosmology-dependent covariance 

Our grafted covariance of the two-point correlation function is es- 
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timated for a fiducial cosmological model, which is given by the 
iV-body simulations. In order not to bias the likelihood function 
of the data (Sect. 14. 2t at points other than that fiducial model, we 
need to account for the fact that the covariance depends on cos- 
mological parameters. We mo del the cosmology -dependence of 
the covariance matrix following Eifl er et al. N2009h . who suggested 
approximative schemes for the mixed term M and the cosmic vari- 
ance term V. Accordingly, for the cosmic-variance term, we as- 
sume a quadratic scaling with the shear correlation function. This 
is true on large scales, where the shear field is close to Gaussian 
and the covariance is indeed proportional to the square of the cor- 
relation function. We calibrate the small-scale Clone covariance in 
the same way, as any differences in the way the non-Gaussian part 
might scale are likely to be small. 

For the mixed term M, we use the fitting formula provided by 
Eifle r et al. I d2009l) . They approximate the variation with S7 m and 
as, leaving the matrix fixed for other parameters. The shot-noise 
term D does not depend on cosmology. The final expression for the 
covariance matrix is 



/(p) 



Di$ij Sj_iv 



+ AV,y(po)( ^ 



0.9/ 



(12) 



where the indices a, v stand for the components '+' and '— ', and 
i, j are the angular-scale indices. Here, p = (fi m , <rg, . . .) denotes 
the cosmological parameter for which the covariance is evaluated, 
with po being the fiducial model of the Clone simulation (Sect. 
13.31 ). The shot-noise term Di and the mixed term M U y,i,i a re es- 
timated using the method of iKilbinger & Schneidej d2004r) . The 
cosmic-variance term Vpu.ij is given in eq. i ll It . The power-law in- 
dices a and j3 depend on the angular scales and the cova riance com- 
ponen t {(j,, v) G {'+',' — '}, and have been obtained in lEifleret alj 
d2009h . We note that the inverse covariance is very sensitive to two 
apparent outliers of a and j3 for the C+--part. To avoid numerical 
issues, we replace these two numbers by the mean of their neigh- 
bouring values. In our case the mixed ter m of the covariance is 
important and cannot be neglected (see also lKilbinger & Schneidei 



l2004IVafaei et alfeoioh . in contrast to recent findings bv lJee et al 
fe012) . 

During the Monte-Carlo sampling (Sect. |4~2t , the covariance 
is updated at each sample point p using eq. d!2t . We make sure 
that each calculation of the covariance resulted in a numerically 
positive-definite matrix, and discard the (rare) sample points for 
which this is not the case. 

In Fig. [4] we show the total covariance C = D + M + V and 
compare it to the Gaussian prediction Cc; = D + M + Vq. Both 
cases are similar on most scales. On small scales the grafted covari- 
ance shows stronger cross-correlations between scales, indicating 
non-Gaussian effects. We find that the additional covariance C m 
due to the shear calibration (see Sect. l3.4t can be neglected, as can 
be seen in Fig.[3]and Sect. 16.21 



3.3.3 Inverse covariance estimator 



It has been shown in lAndersor] d2003l) and lHartlap et all J2007I) that 
the maximum-likelihood (ML) estimator of the inverse covariance 
is biased high. The field-to-field covariance from the Clone is such 
an ML estimate. The bias depends on the number of realisations 
or fields n, and the number of bins p. The ML estimator can be 
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Figure 5. Trace of the covariance and inverse grafted cosmic covariance, 
plotted against the ratio of the number of bins p to the number of real- 
isations n. Shown are the cases for covariance (red squares), the inverse 
(green filled circles), and two cases where the inverse has been multiplied 
with the Anderson-Hartlap factor a, corresponding to the number of bins 
p = 10 (blue open circles), and p = 6 (magenta triangles), causing an over- 
correction in both cases. This shows that we can use the inverse covariance 
estimator without correction (green curve). 



de-biased by multiplic ation with the Ande rson-Hartlap factor a = 
(n - p - 2)/(n - 1) faartlap et al]|2007l) . 

Our final 2PCF covariance, however, is the mixture of an ML 
estimate and analytical expressions. The ML estimator is modu- 
lated with the Gaussian cosmic variance via eq. dl It . to which we 
add the shot-noise and mixed terms, eq. d!2t . We quantify a pos- 
sible bias of the inverse covariance (C++) -1 by varying n for a 
fixed p — 10. For a step index s = 7, corresponding to $ s = 37 
arcmin, Fig. [5] shows that the trace of (C++) -1 does not depend 
on the ratio p/n for our grafted cosmic covariance matrix. Multi- 
plication with a for p = 10 results in an overcorrection, causing 
a strong decrease of tr(C++)~ 1 with p/n. A similar albeit less 
strong decrease is seen when naively taking into account the fact 
that the Clone covariance only contributes to an effective number 
of scales p e g = 6, according to eq. I ll It with s — 7. The three 
curves for the inverse seem to converge for p/n — > 0. Therefore, 
we have reason to be confident that any bias of the unaltered inverse 
of eq. dl2t is small, and hence we do not need to apply the scalar 
correction factor a. The addition of a deterministic component to 
the ML covariance seems to be sufficient to render the estimate of 
the inverse to be unbiased. 



3.3.4 Covariance of derived second-order functions 

Expressions for the covariance of the derived second-order statis- 
tics (eq. [6} are straight-forward to obtain, and can be calculate d 
by integrating the covariance of the 2PCFs l lSchneider etai1l2002h . 
However, the necessary precision for the numerical integration re- 
quires a large number of angular bins for which the 2PCF covari- 
ance has to be calculated, which is very time-consuming. Conse- 
quently, for all derived second-order functions we choose not to 
graft the Clone covariance to the Gaussian covariance, but instead 
only use the Clone to calculate the total covariance of the derived 
functions. To include shot noise, we add to each galaxys' shear an 
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intrinsic ellipticity as a Gaussian random variable with zero mean 
and dispersion o> = 0.38. The latter is calculated as <j\ = £%£*, 
where the sum goes over all CFHTLenS galaxies in our redshift 
range. Therefore, the covariance between the 184 Clone lines of 
sight gives us the total covariance D + M + V. Contrary to the case 
of the 2PCFs (previous section), this covariance stems from a pure 
ML estimate, and therefore the inverse needs to be de-biased by 
the Anderson-Hartlap factor a. With a typical number of angular 
scales of p — 10 to 15 the corresponding a is of order 0.9. We 
show that our cosmological results are independent of the number 
of realisations in Sect. 16.21 Note that for the all derived estimators, 
the cosmology-dependence of the covariance is neglected. 

For upcoming and f uture tomographic su rveys such as KiDS0, 
DES^, HSCfl Euclid^ l lLaureiis et all 1201 ll) or LSST0, a much 
larger suite of simulations will be necessary. The number of re- 
ali sations n has to be substantially larger than the number of bins 
p faartlap et al]|2007t) . For a multi-bin tomographic shear survey, 
p can easily be of the order of several hundreds or more if other 
probes are jointly measured such as galaxy clustering or magnifi- 
cation. This necessitates on the order of a thousand and more inde- 
pendent lines of sight. This number has to be multiplied by many 
if a proper treatment of the cosmology-dependence is to be taken 
into account. Moreover, a simple up-scaling of smaller simulated 
fields to full survey size might not be easy because of the different 
area-scaling of the HSV term. 




a io" 



10" 



Figure 6. The measured shear correlation functions (black squares) and 
£_ (blue circles), combined from all four Wide patches. The error bars cor- 
respond to the total covariance diagonal. Negative values are shown as thin 
points with dotted error bars. The lines are the theoretical prediction using 
the WMAP7 best-fitting cosmology and the non-linear model described in 
Sect. 14. 31 The data points and error bars are listed in Table IbTI 



3.4 Ellipticity calibration corrections 

We ap ply the shear calibration as described in iHevmans et al.l 
d2012l) . which accounts for a potential additive shear bias c and 
multiplicative bias m, 

obs /i i \ true . 

e = (1 + m)e + c. (13) 

The additive bias is found to be consistent with zero for ei. The sec- 
ond ellipticity component £2 shows a signal-to-noise ratio (S/N) 
and size-dependent bias which we subtract for each galaxy. This 
represents a correction which is on average at the level of 2 x 10~ 3 . 
The multiplicative bias m is modelled as a function of the galaxy 
S/ N and size r. It is fit simultaneously in 20 bins of S/N and r, 
see lMiller et all | |2013|) . We use the best-fitting function m(S/N, r) 
and perfor m the global correct ion to the shear 2PCFs, see eqs. (19) 
and (20) of lMiller et afl d2013l) . Accordingly, we calculate the cali- 
bration factor 1 + K as the weighted correlation function of 1 + m, 



, , J2a WiWj(l + m,i)(l + rtij) 
1 + KOff) = =2 ^ -. 

W i W 3 



(14) 



The final calibrated 2PCFs are obtained by dividing and £_ by 
1 + K. The amplitude of 1 + K is around 0.91 on all scales. The 
errors on the correlation function from the fit uncertainty are negli- 
gible compared to our statistical errors. Furthermore, we calculate 
the covariance matrix C m for the correlation function from this un- 
certainty, and show in Sect. l6.2l that the cosmological results remain 
unchanged by adding this term to the analysis. 

Figure|6]shows the combined and corrected 2PCFs, which are 
the weighted averages over the four Wide patches with the number 
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Figure 7. The measured shear correlation functions (top panel) and f _ 
(bottom), for the four Wide patches. The error bars correspond to Poisson 
noise. 



of pairs as weights. Note that the data points are strongly corre- 
lated, in particular on scales larger than about 10 arcmin. Cos- 
mological results using this data will be presented in Sect. [5] The 
correlation signal split up into the contributions from the four Wide 
patches is plotted in Fig. [7] There is no apparent outlier field. The 
scatter is larger than suggested by the Poisson noise on large scales, 
in agreement with the expected cosmic variance. 
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3.5 E- and B-modes 

The aperture-mass dispersion is shown in the upper panel of Fig.[8] 
The B-mode is consistent with zero on all scales. We quantify this 
by performing a null \ 2 test, taking into account the B-mode Pois- 
son covariance C x as measured on the Clone, 

xl = ]T <M X > (Ot) [C- 1 ]^ <Mx) (6,). (15) 

ij 

Since here the covariance is entirely estimated from the Clone line- 
of-sight, the inverse has to be de-biased using the Anderson-Hartlap 
factor. We consider the B-mode over the angular range [5.5; 140] 
arcmin. As discussed before, the lower scale is where the B-mode 
due to leakage is down to a few per cent. The upper limit is given 
by the largest scale accessible to the Clone, which is much smaller 
than the largest CFHTLenS scale: It is 280 arcmin, resulting in an 
upper limit of (Mf p ) of half that scale. The resulting xVdof of 
14.9/15 = 0.99 , corresponding to a non-null B-mode probability 
of 46 per cent. Even if we only take the highest six (positive) data 
points, we find the \ z per degree of freedom (dof) to be x 2 /dof = 
4.12/6 = 0.69, which is less than la significance. The non-zero 
B-mode signal at around 50 - 120 arcmin from F08 is not detected 
here. 

The top-hat shear rms B-mode is consistent with zero on all 
measured scales, as shown in the middle panel of Fig. [8] Note, 
however, that of all second-order functions discussed in this work, 
(|7| 2 ) is the one with the highest correlation between data points. 
The predicted leakage from the B- to the E-mode is smaller than 
the measured E-mode, but becomes comparable to the latter for 
6 > 100 arcmin, where the leakage reaches up to 50 per cent of the 
E-mode. 

The optimized ring statistic for r\ = $ m in/$m ax = 1/50 is 
plotted in the lower panel of Fig. [8] Each data point shows the E- 
and B-mode on the angular range between $ m i n and $ max , the lat- 
ter of which is labelled on the a;-axis. The B-mode is found to be 
consistent with zero, a \ 2 null test yields a 35 per cent probability 
of a non-zero B-mode. 

We first test our calculation of COSEBIs on the CFHTLenS 
Clone with noise, where we measure a B-mode of at most a few 
x 10~ 12 for n ^ 5 and # max *S 250 arcmin. Even though this is a 
few orders of magnitudes larger than the B-mode due to numerical 
errors from the estimation from theory, it is insignificant compared 
to the E-mode signal. When including the largest available scales 
for the Clone however, $ max ~ 280 arcmin, the B-mode increases 
to be of the order of the E-mode. This is true independent of the 
binning or whether noise is added. We presume that this is due to 
insufficient accuracy with which the shear correlation function is 
estimated from the simulation on these very large scales, from only 
a small number of galaxy pairs. Further, for n > 5 a similarly 
large B-mode is found for some cases of (l? m i n , $ max ). Again, the 
accuracy of the simulations is not sufficient to allow for precise nu- 
merical integration over the rapi dly-oscillating filter functions of 
Log-COSEBIs for higher modes ( IBeckenl2012n . We will therefore 
restrict ourselves to n ^ 5 for the subsequent cosmological analy- 
sis. 

The measured COSEBIs modes are shown in Fig. [9] We use 
as smallest scale 1? = 10 arcsec, and two cases of i? max of 100 and 
250 arcmin. In both cases we do not see a significant B-mode. The 
signal-to-noise ratio of the high mode points decreases when the 
angular range is increased: For $ max = 250 arcmin only the first 
two modes are significant. This is not unexpected, since the filter 
functions for i9 max = 250 arcmin sample larger angular scales and 



put less weight on small scales where the signal-to-noise ratio in 
the 2PCFs is larger. 

A further derived second-or der quantity are the shear E-/B- 
mode coiTelation functions £e,b dCrittenden et alj|2002t IPen et alj 
l2002h . which have been used in F08. Whereas they share the in- 
convenience with the top-hat shear rms of a formal upper infinite 
integration limit, they offer no advantage over the latter, and will 
therefore not be used in this work. 

3.6 Conclusion on estimators 

We compared various second-order real-space shear functions, 
starting with the fundamental two-point coiTelation functions £±. 
From the 2PCFs we calculated a number of E-/B-mode separat- 
ing functions. The top-hat shear rms (|7| 2 ) is of limited use for 
cosmological analysis because of the cosmology-dependent E-/B- 
mode leakage. For the aperture-mass dispersion (Aff p ) this leakage 
is confined to small scales, whereas the optimized ring statistic TZe 
and COSEBIs were introduced to avoid any leakage. The drawback 
of the 2PCFs is that they are sensitive to large scales outside the 
survey area and thus m ay contain an undetectable B-mode signal 
jSchneidere7aT]|2010h . COSEBIs capture the E-/B-mode signals 
in an optimal way on a finite angular-scale interval [i9 m i n ; $m ax ]- 
The interpretation of COSEBIs and the matching of modes to an- 
gular scales are not straightforward since the corresponding filter 
functions are strongly oscillating. 

For lensing alone, we obtain cosmological parameter con- 
straints on fl m and as for the different estimators discussed in this 
section. The results and comparisons are presented in Sect. 15. 1 1 

We decided to use the 2PCFs to compute cosmological con- 
straints in combination with the other probes, for the following rea- 
sons. The goal of this paper is to explore the largest scales available 
for lensing in CFHTLenS. This is only possible with a sufficiently 
large signal-to-noise ratio when using the 2PCFs. We note that on 
these lar ge scales our systema tics tests, the star-galaxy shape cor- 
relation dHevmans et al. I I2OI2T) and the E-/B-mode decomposition 
(this work), were not possible. However, since both tests have re- 
vealed no systematics on smaller scales, we are confident that the 
shear signal up to very large scales is not significantly contam- 
inated. Moreover, the implementation of a cosmology-dependent 
covariance is currently only feasible for the 2PCFs. 



4 COSMOLOGY SET-UP 
4.1 Data sets 

We use the following data sets and priors: 

(i) CFHTLenS two-point shear correlation functions and covari- 
ance as described in Sect. [3] We choose the smallest and largest 
angular bins to be 0.9 and 300 arcmin, respectively. This includes 
galaxy pairs between 0.8 and 350 arcmin. 

(ii) Co smic microwave background (CMB) anisotropies: 
WMAP7 dLarson et al]|201 ll ; iKomatsu et alj|201 ll) . The released 
WMAP co dtEB is em ployed to calculate the likelihood, se e also 
iDunklev et alTT2009l) . We use C AME0 dLewis et al.ll200d) to get 
the theoretical predictions of CMB temperature and polarisation 
power- and cross-spectra. 

10 http://lambda.gsfc.nasa.gov 

11 http://camb.info 
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Figure 8. Smoothed second-order functions: Aperture-mass dispersion 
(Af^p) (left panel), shear top-hat rms (|7| 2 ) (middle) and optimized ring 
statistic 7?-e (right), split into the E-mode (black filled squares) and B- 
mode (red open squares). The error bars are the Clone field-to -field rms. 
The dashed line is the theoretical prediction for a WMAP7-cosmology (with 
zero E-/B-mode leakage), the dotted curve shows the Clone lines-of-sight 
mean E-mode signal. For (M 2 p ) and (|7| 2 ) the WMAP7-prediction of the 
leaked B-mode is shown as red dashed curve; the shaded region in the mid- 
dle panel corresponds to the 95 per cent WMAP7 confidence interval of <jg 
(flat ACDM). For the shear top-hat rms, negative points are plotted with 
dashed error bars. 
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Figure 9. COSEBIs with logarithmic filter functions. The left (right) panel 
correspond to a maximum angular scale of 100 arcmin (250 arcmin). The 
filled (open) squares correspond to the CFHTLenS E-mode (B-mode). The 
error bars are the Clone field-to-field rms (rescaled to the CFHTLenS area). 
The dashed line is the theoretical prediction for a WMAP7-cosmology, the 
dotted curve shows the Clone mean COSEBIs. 



(iii) Baryonic acoustic oscillations (BAO): SDSS-III (BOSS). 
We use the ratio D v fr s = 13.67 ± 0.22 of the apparent BAO 
at 2 = 0.57 t o the sound horizon di stance, as a Gaussian random 
variable, from lAnderson et al ] d2012l) . 

(iv) Hubble constant. We add a Gaussian prior for the Hubble 
constant of h = 0.742 ± 0.036 from Cepheids an d nearby type-la 
supernovae distances from HST dRiess et alj|2009l . hereafter R09). 



In contrast to lKilbinger et all d2009h we do not include super- 
novae of type (SNIa) la. BOSS puts a tight constraint on the expan- 
sion history of the Universe, which is in excellent agreement with 
corresponding constraints using t he luminosity distance from the 
most recent compilation of SNIa dConlev et alj|20 1 ll : ISuzuki et all 
l20H) . Both BOSS and SNIa are geometrical probes, and adding 
SNLS to WMAP7+BOSS yields little improvement on cosmolog- 
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ical p arameter constraints with the exception of w dSanchez et alj 
12012) . 

All data sets are treated as independent, neglecting any covari- 
ance between those probes. Experiments which observe the same 
area on the sky are certainly correlated since they probe the same 
cosmological volume. However, this is a second-order effect, like 
CMB lensing, the integrat ed Sachs-Wolfe effect (ISW), or the lens- 
ing of the baryonic peak dVallinotto et~al]|2007l) . Compared to the 
statistical errors of current probes, these correlations can safely be 
ignor ed at present, but have to b e taken into account for future sur- 
veys dGiannantonio et al]|2012h . 



4.2 Sampling the posterior 

To obtain constraints on cosmological parameters, we estimate the 
posterior density n(p\d, M) of a set of parameters p, given the data 
d and a model M. Bayes' theorem links the posterior to the likeli- 
hood L(d\p, M), the prior P{p\M) and the evidence E(d\M), 

L(d\p,M)P(p\M) 



n(p\d, M) 



(16) 



E(d\M) 

To estimate the true, unknown likelihood distr ibution L, a suite of 
JV-body simulati ons would be necessary (e.g. lHarflap et al.ll2009l ; 
IPires et al]|2009h . This is not feasible in a high-dimensional pa- 
rameter space, and for the number of cosmological models probed 
in this work. Instead, to make progress, we use a Gaussian likeli- 
hood function L, despite the fact that neither the shear field nor the 
second-order shear functions are Gaussian random fields. Never- 
theless, this is a reas onable approxim ation, in particular when CMB 
is added to lensing dSato et alfeoioh . 

The construction towards the true likelihood function can be 
informed by further features o f the estimators, for examp le con- 
strained correlation functions dKeitel & Schneider! 1201 ll These 
constraints are equivalent to the fact that the power spectrum is pos- 
itive. We do however not attempt to make use of these constaints. 
The expected deviations are minor compared to the statistical un- 
certainty of the data. The likelihood function is thus given as 

L(d|p,A/)=(2^)- m/2 !C(p,M)r 1/2 

x exp [(d - y( P , M)f C-'ip, M) (d - y(p, M))} , 

(17) 

where y(p, M) denotes the theoretical prediction for the data d for 
a given m-dimensional parameter vector p and model M. 

We sample the posterior with Popul ation Monte Carlo (PMC; 
IWraith et al.ll2009llKilbinger et alj|20ich. using the publicly avail- 
able code COSMO.PM d 12 l dKilbinger et alj|201 lh. PMC i s anadap- 
tive importance-sampling technique dCappe et al.1 120041 120081) in 
which samples p n , n = 1 ... N are created under an importance 
function, or proposal density q. The sample can be used as an es- 
timator of the posterior density ir, if each point is weighted by the 
normalised importance weight 



?(Pn) ' 



ton 



(18) 



The main difficulty for importance sampling is to find a suitable 
importance function. PMC remedies this problem by creating an 
iterative series of functions qt,t — 1 ... T. In each subsequent iter- 
ation, the importance function is a better representation of the pos- 
terior, so the distribution of importance weights gets progressively 



narrower. A measure for this quality of the importance sample is 
the normalised Shannon information criterion, 



Hn = — 53 w n log i 



(19) 



As a stopping criterion for the PMC iterations, we use the related 

perplexity p, 



p = exp(H N )/N, 



(20) 



which lies between and 1, where 1 corresponds to maximum 
agreement between importance function and posterior. 

Most PMC runs reach values of p > 0.7 after 10 or 15 it- 
erations. To obtain a larger final sample, we either perform a last 
importance run with five times the number of points, sampled un- 
der the final importance function, or we combine the PMC samples 
with the five highest values of p. In each iteration we created 10k 
sample points; the final sample therefore has 50k points. 

An estimate E of the Bayesian evidence 



E 



d m pL(d\p,M)P(p\M) 



(21) 



is obtained at no furt her computing cost from a PMC simulation 
lKilbingeretai]d20ld) , 



E 



-T' 



N 



(22) 



4.3 Theoretical models 

We compare the measured second-order shear functions to non- 
linear models of the large-scale structure, with a prediction of the 
density power spe ctrum from the halofit fitting formulae of 
ISmith et alj d2003h . Fo r dark-energy models , we adopt the scheme 
of the ICOSMc(3 co de dRefregier et al.l2008l) . which uses the open- 
CDM fitting formula for a model with wo = —1/3, and inter- 
polates between this case and AC DM for models w i th diff ering 
wo. This scheme was employed in ISchrabback et all d20ld) who 
compa red their non-linear power spectrum with tvlcDonald et alj 
d2006h and found good agreement in the range wq £ [—1.5; —0.5] 
out to k of a few inv erse Mpc in the relevant redshift range. 
Vanderveld et al. I J2012h have shown that for ACDM the halofit ac- 
curacy of 5 to 10 per cent is sufficient for current surveys. From 
hydro-dynamical simulations, baryonic effects have been quanti- 
fied. The results depend on the scenario and specific baryonic pro- 
cesses included in the simulations. The bias in the power spec- 
trum to k of few inverse Mpc is between 10 and 20 per cent. 
The resulting bias on cosmolo gical parameters is smaller than the 
CFHTLenS statistical errors jjing et al.1 120061; iRudd et al.1 120081; 
ISembolonietalJl201lh . 

We also see a good agree ment with the ACDM simulations 
of lHarnois-Deraps et al.1 d2012l) . A more accurate non-linear power 
spectrum on a wider range of cosmo logical parameters could be 
obtained from the Coy ote emulatoi F^I dHeitmann et alj|2010l . l2009l ; 
lLawrence et al. I l2010u . Unfortunately it is limited in wave mode 
(fc < 2.4 Mpc) and, more importantly, by an upper redshift of 
z = 1. Due to the scatter in photometric redshifts, we would have 
to cut at a very low redshift to limit the redshift tail at z > 1. For 
example, with z p h ^ 0.8, the fraction of galaxies at z > 1 is 5 



www . cosmopmc . info 



" www.icosmo.org 

14 http : / /www . lanl . gov/pro ject s/ cosmology/ CosmicEmu 
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per cent, and ignoring these galaxies would bias low the mean red- 
shift by 0.05 which would result in a bias on erg which is larger 
than the uncertaint y in the halo fit prescription. Alternatively, the 
hybrid approach of lEiflerl j201 ll) could be taken, which pastes the 
halofit power spectrum to the Coyote emulator outside the k- 
and z-validity range of the latter. However, this implies multiplying 
one of the spectra by a constant to make the combined power spec- 
trum continuous, We do not deem this sufficiently justified, since 
this multiplicative factor does not stem from a fit to numerical sim- 
ulations and might introduce a bias to erg. 

We individually run PMC for CFHTLenS and WMAP7, re- 
spectively. For the combined posterior results, we importance sam- 
ple the WMAP7 final PMC sample, multiplying each sample point 
with the CFHTLenS posterior probability. 

For weak lensing only, the base parameter vector for the flat 
ACDM model is p = (fi m , erg , fib , n s , h) . It is complemented by 
u>o and fide for dark-energy and non-flat models, respectively. With 
CMB, we add the reionisation optical depth r and the Sunyaev- 
Zel'dovich (SZ) template amplitude Asz to the parameter vector. 
Moreover, we use as the primary normalisation parameter, and 
calculate erg as a derived parameter. We use flat priors through- 
out which, when WMAP7 is added to CFHTLenS, cover the high- 
density regions and the tails of the posterior distribution well. 

For model comparison, we limit the parameter ranges to 
physically well-motivated priors for those parameters which vary 
between models. This is important for any interpretation of the 
Bayesian evidence, since the evidence directly depends on the prior. 
The prior is an inherent part of the model, and we want to compare 
physically well-defined models. 

Thus, we limit the total matter and dark-energy densities fi m 
and fide £ [0; 1], setting a lower physical limit, and creating a 
symmetrical prior for the curvature Qk of [—1; 1], which is bound 
from below by the physical limit of an empty universe. Note that 
by sampling both fi m and fide, the curvature prior is no longer 
uniform but has triangular shape. 

For the model comparison cases we limit wq to [—1; —1/3], 
therefore excluding phantom energy and dark-energy models which 
are non-accelerating at the present time. These priors are the same 
as for the models t hat w ere compared using the Bayesian evidence 
in iKilbinger etaLI < l2010h . The prior ranges for the other parameters 
are fi b G [0;0.1],r G [0.04; 0.2], n s G [0.7; 1.2], 10 9 A|, G 
[1.8; 3.5], h G [0.4; 1.2] and A S z G [0;2]. For the dark-energy 
model runs for parameter estimation, which are not used for model 
comparison, we use a wide prior on wo which runs between —3.5 
and 0.5. 



5 COSMOLOGICAL RESULTS 

The most interesting constraints from 2D weak lensing alone are 
obtained for fi m and erg, which we discuss below for the four cos- 
mologies considered here. Table [2] shows constraints from lens- 
ing alone on the combination erg(fi m /0.27) a , which is the di- 
rection orthogonal to the fi m - erg degeneracy 'banana' To ob- 
tain a, we fit a power law to the log-posterior values using his- 
togr ams with opt imal bin numbers for estimating the posterior den- 
sity JScotdl 19 79h. (Fig. UOt . We also discuss constraints on fiA 
(for cases with free curvature) and wo (for uiCDM models). Table 
[3] shows the combined constraints from CFHTLenS+WMAP7 and 
CFHTLenS+WMAP7+BOSS+R09. The comparison between cos- 
mological models is shown in Tableland described in Sect. 15.41 



5.1 fi m and erg 

Flat ACDM For a flat ACDM universe, the constraints in the 
fim — erg plane (left panel of Fig. [TO) from CFHTLenS are nearly 
orthogonal to the ones for WMAP7. CFHTLenS improves the joint 
constraints for these parameters by a factor of two. Lensing plus 
CMB constrain fi m and erg to better than 5 per cent and 2 per cent, 
respectively. Adding BOSS and R09 decreases the error on fi m to 
3.5 per cent, but does not improve the constraint on erg. 

Flat toCDM If the dark-energy equation-of-state parameter wo is 
kept free, CMB and lensing display the same degeneracy direction 
between fi m and as (left panel of Fig.ll It. Combining both probes 
only partially lifts this degeneracy, the uncertainty on fi m remains 
at the 25 per cent level. This uncertainty decreases to 10 percent 
with the addition of the BOSS BAO distance measure. 

The value of the Hubble constant from both 
CFHTLenS+WMAP7 (h = O.ffit^H) and BOSS+WMAP7 
(h = 0.65iQ'o|) are slightly lower when compared to the R09 
result, h — 0.742 ± 0.036, although it is within the 1-er error 
bar. Since h is degenerate with all other parameters except n s , 
those parameter means change with the inclusion of the R09 prior. 
This causes the relatively large fi m and fib and low erg if R09 
is not added. The joint Hubble constant with all four probes is 



Curved ACDM With curvature left free and no additional priors, 
CMB anisotropies cannot determine fi m anymore, since there is a 
degeneracy between matter density, curvature and the Hubble con- 
stant. Lensing however shows a similar dependency on fi m and erg 
to the flat model case. Therefore, the improvement on fi m from 
CFHTLenS + WMAP7 with respect to WMAP7 alone is an order 
of magnitude, to yield a 8 per cent error. The joint error on erg is 
3.5 per cent. 

Curved wCVM The fi m - erg degeneracy holds nearly the same 
as in the previous cases of models with fewer parameters, as dis- 
played in the left panel of Fig. [12] The value of crg(fi m /0.25) a 
is slightly increased but well within the error bars. The joint 
CFHTLenS+WMAP7 results on fi m and erg are similar to the flat 
toCDM case. 

The BOSS+R09+WMAP7 results indicate a slightly smaller 
fi m and larger erg. The joint CFHTLenS+WMAP7+BOSS+R09 
allowed region is therefore on the upper end of the 
CFHTLenS+WMAP7 banana. The reason for this is, as in 
the flat wCDM case, the degeneracy of fi m and erg with the Hubble 
constant. WMAP7 alone prefers a low value, h = 0.51q'i3, which 
increases to h = 0.73 ± 0.04 when BOSS+R09 is added. As a 
consequence, fi m decreases and erg increases. On the other hand, 
adding CFHTLenS to WMAP7 leaves the Hubble constant at the 
relatively low value of h = 0.60t^;^. 

5.2 Dark energy 

For the following results on the dark-energy equation-of-state pa- 
rameter w, we use the flat prior [—3.5; 0.5]. 

Flat toCDM 2D weak gravitational lensing alone is not able to 
tightly constrain dark energy, in contrast with 3D tomographic 
weak lensing. The 68 per cent confidence limits for wo (flat 
mCDM) are of order unity, Wo = — *|. In combination with 
WMAP7 only, these errors decrease by a factor of four, and wo gets 
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Table 2. Constraints from CFHTLenS orthogonal to the f2 m -CTg degeneracy direction, using the 2PCF. The errors are 68% confidence intervals. The four 
columns correspond to the four different models. 



Parameter 


flat ACDM 


flat tuCDM 


curved ACDM 


curved tuCDM 


<7 8 (fW0.27)« 
a 


0.79 ± 0.03 
0.59 ±0.02 


7O+0.07 
u - ' y -0.06 

0.59 ± 0.03 


0.80±°;« 
0.61 ±0.02 


0.61 ±0.03 



Table 3. Cosmological parameter results with 68% confidence intervals. The first line for each parameters shows CFHTLenS+WMAP7, the second line is 
CFHTLenS+WMAP7+BOSS+R09. The star (* ) indicates a deduced parameter. The four columns correspond to the four different models. 



Parameter 


flat ACDM 


flat u>CDM 


curved ACDM 


curved miCDM 




274+ - 013 
283 +0 1 ' 01 ,„ 


n o 9 c+0.082 

u.ozo_ 076 
u.^oi _Q.023 


n 97 ,-+0.023 
u - z '°-0.021 

286+?'?, 1 , 1 


n 077+0. 098 
' -0.079 

271+ ' . 2 , 8 

^-0.025 


<X 8 * 


n 81^+0-016 
0-810_o .014 

814+2-S5 

u - ol *-0.014 


77+ - 11 
u -' '-0.07 

809+S ^? 

u.oua_Q Q 35 


n 81 ej+0.030 
U-»10_o .025 
804 + 018 


n 7ic+0.090 
u. / rJ_o.o 7 Q 

U ' OZU -0.039 


Wo 


-1 


-0.86l ;| 
-0 99+2" 


-1 


n 79 +0.20 
u - ,z -0.24 

-1 10+°-" 

Xi±u -0.16 




i-n m 


i - n m 


0.726i° if 
0.7186I™? 9 


0.628±g : g« 

n 7OCT+0.028 
U. /OO_ 032 


n K * 








-0.0003 ± 0.0086 
-o nn47 +0 ' 0045 

u - uu4 ' -0.0047 


-0 00 ,+ alm 
u.uuo_ 012 

-0 0063+ 0064 
u.uuoj_ 0045 


h 


702+ ' 014 
u -' uz -0.013 

0.693 ±0.010 


0.66lg; 07 
0.69ll ;gi 


n 7n o+0.037 
u. (ua_ 033 

0.683 ±0.014 


605+ - 082 
u.ouo_ 062 

7n2+ - 032 
u. /uz_ 030 


n h 


0.0456 ± 0.0012 
0.0465 ± 0.0010 


U.UO4_ 013 

n471+ - 0046 
u.U4/i_ 0042 


n457+ - 0045 

U.U40C_q 0Q41 

n n489+°' 0020 
U.0482_ 0019 


064+ - 018 
u.uo4_ 014 

0457+ - 0047 




0.966 ±0.013 
0.961 ±0.012 


0.966 ±0.014 
0.9591™ 1 ! 


qfi5+ 013 
u.aoo_ 014 

0.966 ±0.013 


Q70+ - 014 
u.y/u_ 013 

0.964I™ 13 


T 


o.o89i° :° is 

0-083±g : glS 


0.088t ™ 
0.084l°J 15 3 


O.OSSlHi! 
0.088±Hi 6 4 


0.0881™" 
0.0871°:™ 


10 9 A| 


2.44il°;°t2 

9 4 = 7+0.088 
^■40 1 -0.081 


n .oq+0.095 

2 465+0- 097 
Z -^ DO -0.089 


9 499 + 0.0H5 
z '^ zz -0.088 


9 q Q c+0.093 
Z -^ ZO -0.089 



Adding BOSS+R09 decreases the error by another factor of two to 
around 0.005. The combined constraints are thus consistent with a 
flat universe within 5 x 10 -3 . For a tuCDM model, this uncertainty 
is of the same order. 



5.4 Model comparison 

In Table[4] the evidence E and the logarithms of the evidence ratios, 
lnBoi = IxiEq/Ei between the baseline flat A model and the 
other three models are shown. Here, Eo is the evidence for flat 
ACDM and E\ the evidence for one of the three models curved 
ACDM, flat iuCDM, and curved u>CDM, respectively. Box is called 
the Bayes factor between model '0' and model T. 

An empirical scale to interpr et such eviden ce values was sug- 
gested by |jeffrevsl ( fl96lh . see also lTrottal l llooih . Accordingly, two 
models are not distinguishable when | In B01 1 < 1- If the log-Bayes 
factor is between 1 and 2.5, the evidence is called weak. Moderate 



constrained to about 30 per cent. The CFHTLenS+WMAP7+BOSS 
constraints on dark energy are wo = — 0.781Q 1 ] 19 . We discuss this 
deviation from ACDM in Sect. [7] Adding the R09 prior on H 
does not reduce the error but shifts the mean to the ACDM value, 

w = -0.99lg;^. 

Curved miCDM The case of dark energy is similar in the 
curved case. CFHTLenS alone results in roo = — 1.2l°'|. 
Adding WMAP7 reduced this uncertainty to 30 per cent. 
CFHTLenS+WMAP7+BOSS yield w = -0.8ltg| 9 . Adding the 
R09 prior on Ho, we find the ACDM-consistent value of u>o = 
-1 10+ 015 

1 - iu -0.16- 



5.3 Curvature 

CFHTLenS helps to improve the constraint on the curvature den- 
sity Q.K- For ACDM, the uncertainty decreases by a factor of 10 
from around 0.1 (WMAP7 alone) to 0.01 (CFHTLenS+WMAP7). 
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Table 4. Bayesian evidence E and Bayes' factor with respect to ACDM, for four different cosmological models. The columns for wq and Q.^ show the prior 
range for those two parameters. The data is CFHTLenS+WMAP7 (fourth column) and CFHTLenS+WMAP7+BOSS+R09 (last column), respectively. 



Name 


w 




lnBoi 


lnBoi 








CFHTLenS+WMAP7 


CFHTLenS+WMAP7 










+BOSS+R09 


ACDM 


-1 











curved ACDM 


-1 


[0;2] 


-3.84 


-4.0 


fiat wCDM 


[-1; -1/3] 





0.42 


0.58 


curved mjCDM 


[-1; -1/3] 


[0;2] 


-3.19 


-4.8 



Table 5. Constraints from CFHTLenS orthogonal to the f2 m -<rg degeneracy 
direction. The main result from the 2PCF (first row) are compared to other 
estimators. 



Data 



(7 8 (n m /0.27) c 



2PCF 

<M? P > 

(|7| 2 ) (ignoring offset) 
(I7I 2 ) (constant offset) 

n E 

COSEBIs (i9 max = 100') 
COSEBIs (i? max = 250') 



0.59 ±0.02 
0.70 ±0.02 
0.60 ±0.03 
0.58 ±0.03 
0.56 ±0.02 
0.60 ±0.02 
0.64 ± 0.03 



0.79 ± 0.03 
0.79 ± 0.06 



0.78 
0.80 
0.80 



+0.04 
-0.05 

+0.03 
0.04 

+0.03 
0.04 



79+ ' 04 



0.77 



-0.04 
-0.05 



2PCF, constant covariance 


0.60 


± 


0.03 


n 78 +0.03 


2PCF (-& > 17') 


0.65 


± 


0.02 


0.78 ±0.04 


2PCF (1? > 53') 


0.65 


± 


0.03 


n 70+0-07 

Ui ' 9 -0.06 



evidence is assumed for 2.5 < ] In B01 1 < 5, and strong for values 
larger than that. 

We compute the evidence for the two cases of probes 
CFHTLenS +WM AP7 and CFHTLenS+WMAP7+BOSS+R09. Al- 
though the evidence for the flat dark energy model is slightly larger 
than the one for the cosmological constant model, the models are 
indistinguishable: Their respective evidence values, or posterior 
odds, are within a factor of two. The evidence against curved mod- 
els is moderate, with log-Bayes factor ratios between 3.2 and 4.8, 
or posterior odd ratios between 25 and 130 in favor of flat ACDM. 
The significance increases when adding BOSS and R09, but stays 
in the moderate range. We remind the reader that the dark-energy 
models considered here have the flat prior [— 1 ; — 1/3] for to, which 
corresponds to an accelerating non-phantom dark energy compo- 
nent. 



6.1 Derived second-order functions 

As expected, the constraints from the derived second-order esti- 
mators are less tight than from the 2PCFs, since they always in- 
volve information loss. Moreover, we use a smaller range of angu- 
lar scales, cutting off both on the lower and higher end, as discussed 
before. All estimators give consistent results. 

Aperture-mass dispersion and top-hat shear rms give very sim- 
ilar constraints compared to the 2PCFs. The position and slope of 
the banana are nearly identical, although the width is larger by a 
factor two (see Table[2). For (J7I 2 ), we analyse two approaches of 
dealing with the finite survey-size E-/B-mode leakage: 

(i) Ignoring the leakage. We fit theoretical models of the top-hat 
shear rms (eqs.[7][A3} directly to the measured E-mode data points 
(|7| 2 }(#i). Since power is lost due to the leakage, we expect as fiS 
to be biased low. 

(ii) We add a constant offset of 5.3 x 1CP 7 to the measured E- 
mode points. This corresponds to the theoretical leakage for the 
WMAP7 best-fitting ACDM model with a 8 = 0.8. On scales 
8 < 5 arcmin, the assumption of a constant offset is clearly wrong; 
however, the constant is two orders of magnitudes smaller than the 
measured signal and does not influence the result much. 

The difference between both cases is about half of the statisti- 
cal uncertainty (Table[2}. More sophisticated ways to deal with this 
leakage, e.g. going beyond a constant offset, or marginalising over 
a parametrized offset, are expected to yield similar results. Since 
they all have the disadvantage of depending on prior information 
about a theoretical model which might bias the result towards that 
model, we do not consider this second-order estimator further. 

The function TZe on scales between $ max = 7.5 and 140 ar- 
cmin, for 77 = tfmm/tfmax = 1/50 (implying i? m i n = 9 arcsec 
. . . 2.8 arcmin) yields consistent results to the 2PCFs. The results 
for COSEBIs are consistent with the 2PCFs for both cases of i? max . 
For $ max = 250 arcmin we find that the modes which are con- 
sistent with zero do contain information about cosmology. When 
only using the first two modes, we obtain <78(fi m /0.27) 0,7 = 
0.78^0 o7> corresponding to a larger uncertainty of a 30 per cent. 



6 COMPARISON OF WEAK-LENSING STATISTICS, 
SYSTEMATICS AND CONSISTENCY TESTS 

In this section we obtain cosmological constraints from the derived 
second-order estimators which were discussed in Sect. 12. 21 The fol- 
lowing tests are all performed under a flat ACDM model. The re- 
sults are listed in Table [5] 



6.2 Robustness and consistency tests 

In this section, we test the robustness of our results, by consider- 
ing various potential systematic effects, and by varying the angular 
scales and estimators. 

Shear calibration covariance We add the shear calibration C m 
(Sect. [3~4t covariance to the total shear covariance. The correlation 
coefficient of C m between angular bins is nearly unity, implying 
that the shear calibration varies very little with angular separation. 
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Figure 10. Marginalised posterior density contours (68.3%, 95.5%, 99.7%) 
for CFHTLenS (blue contours), WMAP7 (green), CFHTLenS+WMAP7 
(red) and CFHTLenS+WMAP7+BOSS+R09 (black). The model is flat 
ACDM (left panel) and curved ACDM (middle and right panel), respec- 
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Figure 11. Marginalised posterior density contours (68.3%, 95.5%, 99.7%) 
for CFHTLenS (blue contours), WMAP7 (green), CFHTLenS+WMAP7 
(magenta) and CFHTLenS+WMAP7+BOSS+R09 (black). The model is 
flat toCDM. 



Since the magnitude of the covariance is much smaller than the 
statistical uncertainties, the cosmological results are virtually un- 
changed. 

Large scales only. The largest ratio of signal-to-noise for cosmic 
shear is on small, non-linear scales. Unfortunately, those scales are 
the most difficult to model, because of uncertainties in the dark- 
matter clustering, and baryonic effects on the total power spectrum. 
To obtain more robust cosmological constraints, we exclude small 
scales from the 2PCFs in two cases, as follows. First, we use the 
cut-off # c = 17 arc minutes. At this scale, the non-linear halofit 
prediction of i s within 5 per cent of th e linear model. Baryonic 
effects, following ISemboloni et al. l|), are reduced to sub per 
cent level. The component £_ , being more sensitive to small scales, 
is still highly non-linear at this scale. However, since most of the 
constraining power is contained in £+, the resulting cosmological 
constraints will not be very sensitive to non-linearities. Neverthe- 
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Figure 12. Marginalised posterior density contours (68.3%, 95.5%, 99.7%) 
for CFHTLenS (blue contours), WMAP7 (green), CFHTLenS+WMAP7 
(magenta) and CFHTLenS+WMAP7+BOSS+R09 (black). The model is 
curved uiCDM. 



Figure 13. Marginalised posterior density contours (68.3%, 95.5%, 
99.7%) for WMAP7 (green), WMAP7+CFHTLenS (magenta), 
WMAP7+BOSS+R09 (orange) and WMAP7+CFHTLenS+BOSS+R09 
(black). The model is flat ACDM (upper panel) and curved ACDM (lower 
panel), respectively. 



less, we use a second, more conservative, cut-off of $ c = 53 ar- 
cmin, where the non-linear models of is within a factor of two of 
the linear one. On these scales, is affected b y baryonic physics 
by less than 5 per cent dSemboloni et aljfeoi ll) . In both cases, we 
obtain a mean parameter value for crg(n m /0.27) ' which is con- 
sistent with the result from all angular scales down to an arcmin. In 
comparison, the error bars on this combined parameter are larger 
by 30 per cent for # c = 17 arcmin, and 100 per cent for # c = 53 
arcmin (see Table . 

Reduced shear Since the weak-lensing observable is not the shear 
7, but the reduced shear g = — k), the relation between 

the shear correlation function and the convergenc e power spectrum 
ignor es higher-order terms (see for an overview iKrause & Hiratal 
|2010|) . The full calculation of only the third-order terms, involving 
the convergence bispectrum, is very time-consuming and unfeasi- 



ble for Monte-Carlo sampling, requiring the calculation of tens of 

thousands of different models. 

Instead, we explore the fitting formulae from iKilbingerl ( T20 1 Ot) 
as a good approximation of reduced-shear effects. For a WMAP7 
ACDM cosmology, the ratio between the 2PCFs with and without 
taking into account reduced shear is 1 per cent for and 4 per cent 
for at the smallest scale considered, 1? = 0.8 arcmin. Since the 
fitting formulae are valid within a small range around the WMAP7 
cosmology, we use them for the combined Lensing+CMB parame- 
ter constraints. The changes in fi m and erg for a ACDM model are 
less than a per cent. 

Number of simulated lines of sight Following iHuff et aT] 4201 lh . 
we examine the influence of the number of simulated lines of sight 
on the parameter constraints. We calculate the covariance of (Mf p ) 
from 110 instead of 184 lines of sight (Sect. [333). Using the cor- 
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responding Anderson-Hartlap factor a, we find identical results as 
before and conclude that the number of simulations is easily suffi- 
cient for this work. 

Non-Gaussian covariance fitting formulae 

We replace the non-Gaussian covariance of calculated from 
the Clone with the fitting formulae from ISemboloni et al.l d2007l . 
S07) and I Sato et al.l d201 lL Sll), respectively. These works pro- 
vide prescriptions of the non-Gaussian covariance by rescaling the 
Gaussian cosmic variance term, by fitting to TV-body simulations. 
Since no recipe for is given, we use the Gaussian covariance 
for £_ and for the cross-covariance between and . We find 
for o- 8 (fi m /0.27) Q the results 0.79t°m for the S07 case, and 
0.781q'q4 for Sll, recovering the mean value with slightly larger 
error bars. This shows that our results are not sensitive to the choice 
of the non-Gaussian covariance. 

Cosmology-independent covariance 

We compare the two cases of cosmology-dependent covariance 
(Sect. [3~T2l>, and a constant covariance, fixed to the fiducial model. 
Contrary to ljee et all d2012l) . we find only small differences in the 
cosmological results. The main effect is a slight increase in the error 
bars for the constant covariance, see Table [2] In particular, around 
the region of the fiducial model, the results are basically the same, 
and therefore, the joint constraints with other probes are virtually 
unaffected by the choice of the covariance. 



7 DISCUSSION AND CONCLUSION 

In this paper we present measurements of various second-order 
shear correlations from weak gravitational lensing by CFHTLenS, 
the Canada-France Hawaii Telescope Lensing Survey. Using a sin- 
gle redshift bin, 0.2 < z p < 1.3, we obtain cosmological con- 
straints on the matter density, fi m , and the power-spectrum ampli- 
tude, as. Adding WMAP7, BOSS and R09 data, we obtain parame- 
ter constraints for flat and curved ACDM and dark-energy models, 
and calculate the Bayesian evidence to compare the probability for 
each model given the data. 

Second-order shear functions Along with the two-point correla- 
tion functions and , which are the fundamental shear observ- 
ables, we consider various derived second-order functions, which 
are able to separate the shear correlation into its E- and B-mode. 
The resulting B-mode is consistent with zero on all scales. The ex- 
cess in the E- and B-mode signal that was seen in the CFHTLS- 
T0003 data (F08) between 50 and 130 arcmin is no longer present. 
This excess was most likely due to systematics in the earlier data, 
and the removal of this feature has to be seen as a success of the 
CFHTLenS analysis. In particular, hints for deviatio ns from Gen- 
eral Relativity using the F08 data fchao et alj|2010h are not con- 
firmed with CFHTLenS dSimpson et alj|2013l) . 



Cosmological parameters The parameter combination which 
2D weak lensing can constrain best, is agfl^j with a ~ 
0.6. CFHTLenS alone, with the two-point correlation functions 
(2PCFs), constrains the combination cr 8 (fi m /0.27) ' 6 to 0.787 ± 
0.032. To facilitate a comparison with F08, we write our con- 
straints as a 8 (fi m /0.25) a = 0.82 ± 0.03 (2PCFs), 0.84±g;g| 
(top-hat shear rms) and 0.83 ± 0.06 (aperture-mass dispersion). 
The exponent a is around 0.6 in all three cases. F08 obtained re- 
sults with uncertainties between 0.042 and 0.049, although only 



smoothed second-order quantities were used. The function clos- 
est to the 2PCFs is £ E , for which F08 found cr 8 (fi m /0.25) a46 = 
0.784 ± 0.049, corresponding to an uncertainty of 6.2 per cent, 
compared to 4 per cent in this work. This increase in precision of 
50 per cent is consistent with the naive expectation of CFHTLenS 
to yield smaller error bars by a factor of the square root of the area 
between CFHTLenS and T0003, which is ^129/57 = 1.4. 

Our uncertainty for the top-hat shear rms is about the same 
as in F08, whereas the aperture-mass dispersion is slightly more 
poorly constrained in this work. Tighter constraints of 5.5 per cent 
were found by F08 from (Mf p ), which was used for their combined 
CFHTLS-T0003 + WMAP3 results. The fact that the increase in 
precision for those smoothed quantities is smaller than expected 
is most due to the limited parameter range of F08, who used tight 
priors on the Hubble parameter, h € [0.6; 0.8] and fixed fib and n s . 
In particular the narrow range of h resulted in tighter constraints on 
fi m and as in F08. 

ISchrabback etafl d2010l) obtained a 8 (fi m /0.3) 06 = 0.68 ± 
0.11 from a 2D weak-lensing analysis of the COSMOS data. The 
relatively large error bars are dominated by cosmic variance from 
the very small survey area of 1.64 square degrees, despite the great 
depth of the survey. A 3D lensing analysis with a large number 
of redshift bins (5 narrow bins up to z = 4 and one broad bin 
with mean redshift of 1.5), decreases the error bar by 20 per cent, 
(j 8 (Sl m /0.3) - 5 = 0.79 ± 0.09. This uncertainty is still larger by a 
factor of 3 than our CFHTLenS 2D constraints, CT 8 (fi m /0.3) 0,6 = 
0.74 ±0.03. 

On 168 square degrees of SDSsF^I data in t he Stripe 8 2 equa- 
torial region, out to a median redshift of 0.52, iHuff etal.1 d201ll) 
recently obtained o- 8 (fi m /0.264) ' 67 = 0.65±g;l§ using a combi- 
nation of COSEBIs with ($mi n , $max) = (1.3, 97.5) arcmin and 
an additional data point of (38 arcmin) . With CFHTLenS we get 
(j 8 (fi m /0.264) fi = 0.80 ± 0.03 for the 2PCFs and 0.801^ for 
COSEBIs with i? max = 100 arcmin. 

Recent results from a 2D analysis of the Deep Lens Surve4^l 
(DLS) yielded the very tight con straints fi m — 0.26 ± 0.05 and 
a 8 = 0.87 ± 0.07 Wet alj2012h . Compared to CFHTLenS, DLS 
has a greater depth with mean redshift of 1.1 and 17 galaxies per 
square arcmin, but on the other hand covers with 20 square degrees 
a smaller area. The parameter space sampled by DLS is similar to 
F08, with a tight prior on the Hubble constant and fixed baryon 
density fib and spectral index n s . 

For all models of dark-energy and curvature considered 
here, the agreement of a 8 and fi m from CFHTLenS with 
WMAP7 is very good. This remains true when BOSS data on 
the BAO peak is added. However, we find values of wo for 
CFHTLenS+WMAP+BOSS which are significantly smaller than 
— 1, both for flat and curved uiCDM models. The reason for this is 
the near-degeneracy of the dark-energy parameter with the Hubble 
constant. The latter takes the rather low value of around 0.65 ± 0.1, 
which results in a high value of wo . Adding the R09 result increases 
h and thus also increases wo, yielding values which are consistent 
with ACDM. 

For the flat ACDM model, adding CFHTLenS to WMAP7 
strongly helps reducing error bars on fi m and <r 8 . The improvement 
is larger than in the case where BOSS+R09 is joined with WMAP7, 
in particular for as- The curved ACDM case sees a slightly 
different dependence on Qk between WMAP7+CFHTLenS and 



15 Sloan Digital Sky Survey; www .sdss.org 
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WMAP7+BOSS+R09, resulting in tight constraints when all four 
probes are combined. Both cases are emphazised in Fig. [13] In the 
mCDM case, both CFHTLenS and BOSS cannot improve signif- 
icantly the dark-energy constraint with respect to WMAP7. Only 
the addition of R09, thereby lifting the wo-h degeneracy decreases 
the error on the dark-energy parameter. 

Our results a re in very good agreement with the measure- 
ment presented in iHudson & TurnbulU d2012l) . who find Sl m = 
0.259 ± 0.045 and erg = 0.748 ± 0.035 for a flat ACDM model. 
This method uses low- and high-z peculiar velocity data only and 
is therefore complementary and independ e nt of o ur results. 

Recent constraints by iMantz etal I d2010h from the X-ray 
ROSAT All-Sky Survey using the cluster mass function for a flat 
wCDM universe are Cl m = 0.23 ± 0.04, erg = 0.82 ± 0.05, 
and wo — —1.01 ± 0.20, in agreement with the results pre- 
sented here. Their relatively low fl m is consistent with our result of 
fi m = 0.29±0.02. When adding CMB (WMAP5), S NIa, BAO and 
the clu ster gas fraction to the cluster mass function, iMantz et alj 
j2010l) get fi m = 0.27 ± 0.02. Fro m the optical SDSS maxBCG 
cluster catalogue. lRozo et alJd201ol) obtain for a flat ACDM model 
cr 8 (fi m /0.25) a41 = 0.83 ± 0.03. In combination with WMAP5, 
they get fi m = 0.265 ± 0.016 and erg = 0.807 ± 0.020, which is 
again consistent with this work. 

Model comparison Using the Bayesian evidence, we computed 
the posterior odds for various cosmological models. Starting from 
the basic ACDM model, we tested extensions of this model which 
included curvature Qk and the dark-energy equation of state pa- 
rameter w. We find no evidence against the standard flat ACDM 
model. 

The constraints for the larger models with free curvature are 
consistent with Qk = 0. It is therefore not surprising that those 
more general models are not favoured over models with fixed flat 
geometry. The larger parameter space from the additional degree 
of freedom implies a lower predictive capability of those extended 
models. A good model should not only predict (a priori) the cor- 
rect parameter range where the result is to be found (a posteriori), 
but also make a specific and accurate prediction, in other words, it 
should have a narrow prior range compared to the posterior. A lack 
of predictive capability is penalised by the Bayesian evidence. 

In contrast to the two non-flat models (curved ACDM and 
curved zuCDM), the flat tuCDM universe is indistinguishable from 
a flat model with cosmological constant. This can be understood by 
looking at the respective additional parameter constraints beyond 
ACDM, that is, fi^c for the curved and wq for wCDM. Compared 
to the corresponding prior, the allowed posterior range for fide is 
a lot smaller than the one for wq since the latter parameter is less 
tightly constrained. Therefore, the curved models are less predic- 
tive, corresponding to a lower evidence. Both the very tight con- 
straints on f2_R-, with error of about 0.005, as well as the moderate 
Bayesian evidence in favour of a flat model strengthen the emerg- 
ing picture that we live indeed in a Universe with zero curvature. 
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APPENDIX A: FILTER FUNCTIONS 

We give expressions of the filter functions F + ,F~, needed to cal- 
culate the derived second-order shear observables from the shear 
correlation functions (eq.[6j. See Table Q]for the relation between 
F± and the following functions and the integration ranges. 

Aperture-mass dispersio n The filter functions for the aperture- 
mass dispersion, defined in ISchneideret"ai1 l l2002l) , are for x < 2 



T+(x) = 



6(2 - 15a; 2 



arcsin(a;/2) 



IOOtt 



x (120 + 2320a; 2 - 754a; 4 + 132a; 6 - 9a; 8 



357T 



2 \ 7/2 



(Al) 



The functions have finite support, and are set to zero for x > 2. 
The Fourier-space filter function for the aperture-mass dispersion 
(eq.|7J is 



U(£) = U e {l) 



24J 4 (6^) 



(A2) 



Top-hat shear rms For the top-hat shear rms, the real-space filter 
functions are 

S+(x) =- |4arccos(a;/2) - aV4-a; 2 J H(2 - x) ; 

S-(x) =^-X 

[a;\/4-a; 2 (6 - a; 2 ) - 8(3 - a; 2 ) arcsin(a;/2)], (A3) 

where H is the Heaviside step function. Thus, S+ has support [0; 2] 
wheras S- has infinite support. The Fourier transform of S+ is 

2Ji(0£) 



The function F_ can be obtained by an integral over F+ , which 
follows from the relation eq. [7] see ISchneider et all d2002l) . Apart 
from these conditions, the functions F± can be freely chosen. 

For the optimized ring statistics, the filter functions corre- 
sponding to 7£e,b are linear combinations of Chebyshev polyno- 
mials of the second kind, 



T+(0) =T+(x = 



20 - i? n 



U n (x) = 



JV-l 

: ^2 a„U„(x); 

n=0 

sin[(n + 1) arccos x] 



(A6) 



(A7) 



sin(arccos x) 

The coefficients a„ can be chosen freely, under the condition that 
1Ze,b are p ure E- and B-mode com ponents, respectively. We take 
the a n from lFu & Kilbinge3 d2010h . which minimized the fi m -<78 
la-area using the CFHTLS-T0003 survey setting, and for fixed 
77 = $min/#max = 1/50. For a fixed 77, 1Ze,b depends on only 
one angular scale 6, which we take to be $ max . 

COSEBIs The COSEBIs filter functions we use here are polyno- 
mials in the logarithm of the angular scale 6, 



3=0 



,9 . 
u mil: 



Nn n ( 2 _ r " 

j=i 



(A8) 



U{i) = Ue(£) = 



(A4) 



The p olynomials t l ° g n (z) have been constructed in lSchneider et al.l 
using eq. dA5l as an orthonormal and complete set of func- 
tions. The coefficients c„j are fixed by integral conditions that as- 
sure the E-/B-mode decomposition of the 2PCF on a finite angular 
integral. They are given by a lin ear system of equations, which is 
given in ISchneider etal] J201C|) . To solve this system, a very high 
numeri cal accuracy is needed. W e use the M ATHEMATICA program 
given in Schneider et al. to obtain the coefficients for a given 

i? m in and i9 ma x, and store the zeros r„i, for which a lower accuracy 
is sufficie nt. The function F- is then calculated using eqs. (38) and 
(39) from lSchneider et al.ld2010h . 

Both for COSEBIs and for FK10, no closed expressions for 
the Hankel transforms of T± have been found (yet); neither for 
the Fourier-space counterparts, apart form the linear COSEBIs 
dAsgarietall2012l) . To obtain the theoretical predictions for these 
functions, we first calculate £± via eq. l[4j, and use eq. l[6j to obtain 
the COSEBIs prediction. 

To calculate the numerical integration over the correlation 
function (eq. [6]( with high enough precision, we split up the inter- 
val [0; Zmax] into 10 pieces, and perform a Romberg-integration on 
each piece with relative precision of 10 -6 . The resulting numerical 
B-mode is smaller than 10 -15 for modes n ^ 10, which is about 
three orders of magnitudes smaller than the predicted E-mode. 



Optimized ring statistic 

To obtain an E-/B-mode decomposition of the 2PCF on a fi- 
nite angular range [i? m i n ; i? m ax] via the sum in eq. [6] two in- 
tegral conditions for the filt er function F+ need to be fulfilled 
dSchneider & Kilbinger 200^): 



fl) m „ 

/ dtftf.F + (tf)=/ d00 3 .F+(0) = O. 



(A5) 



APPENDIX B: CFHTLENS SECOND-ORDER 
WEAK-LENSING DATA 

The measured 2PCFs as shown in Fig. [6] We list the data points 
and the total error in Table IB II The full covariance is available on 
request or via the web page http : // cfhtlens . org. The de- 
rived second-order E- and B-mode functions are listed in Tables lB2l 
to|B5] 
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Table Bl. The CFHTLenS two-point correlation functions (2PCFs) and 
£_ , for different angular scales #, see Sect. 13.21 The values of cr indicate 
the error from the total covariance diagonal (Sect. 13.3.1) . The covariance is 
the one at the fiducial Clone cosmology. 



Table B3. The CFHTLenS optimized ring statistic. The E-mode, 7?-e> 
and B-mode 7^.B' are given for different angular smoothing scales 9, see 
Sect. |3.5| The values of cr indicate the error from the scaled Clone covari- 
ance diagonal (Sect. lXX4V 
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Table B4. The COSEBIs for i9 min = 10 arcsec and i9 max = 100 arcmin. 
The E-mode, E n , and B-mode B n , are given for the first five modes n, see 
Sect. 13.51 The values of cr indicate the error from the scaled Clone covari- 
ance diagonal (Sect. l33~4Y 



Table B2. The CFHTLenS aperture-mass dispersion. The E-mode, (Aff p ), 
and B-mode (M x ), are given for different angular smoothing scales 9, see 
Sect. 13.51 The values of cr indicate the error from the scaled Clone covari- 
ance diagonal (Sect. [33!4l . Note that for cosmological results, we do not 
use scales below 5 arcmin. 
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Table B5. The COSEBIs for tf min = 10 arcsec and i? max = 250 arcmin. 
See Table|B4]for details. 
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